Naval  Research  Laboratory 


Washington,  DC  20375-5320 


NRL/MR/5521-95-7744 


Novel  Techniques  for  the 
Analysis  of  Wireless  Integrated 
Voice/Data  Networks 


Jeffrey  E.  Wieselthier 
Craig  M.  Barnhart 

Communication  Systems  Branch 
Information  Technology  Division 


Anthony  Ephremides 

Kaman  Sciences  Corporation 
Alexandria,  Virginia 

and 

University  of  Maryland 
College  Park,  Maryland 


DTIC 

^ELEGTE|% 
^AUG  0  2  1995  1  I 


July  24, 1995 


19950801  014 


Approved  for  public  release;  distribution  unlimited. 


REPORT  DOCUMENTATION  PAGE 

Fom  Approved 

OMBNo.  0704-0188 

Public  reporting  burden  for  this  ooliection  of  titormalion  is  aetimated  to  average  1  hour  per  response.  Including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  correlating  and  reviewing  the  collectbn  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
coMeaion  of  Irrformatlon,  Indoding  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Servicas.  Directorate  for  Information  Operations  and  Reports.  1215  Jefferson 

Davis  Highway.  Suite  1204,  Arlington,  VA  22202-4302.  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188)  Washington,  DC  206503. 

1.  AGENCY  USE  ONLY  (Leave  Bank) 

2.  REPORT  DATE 

July  24. 1995 

3.  REPORT  TYPE  AND  DATES  COVERED 

Interim  Report  6/93-5/95 

4.  TFTLE  AND  SUBTITLE 

Novel  Techniques  for  the  Analysis  of  Wireless  Integrated 

Voice/Data  Networks 

5.  FUNDING  NUMBERS 

PE-61153N 

PR  -  RR015-09-41 

WU  -  DN159-036 

6.  AUTHOR(S) 

Jeffrey  E.  Wieselthier,  Craig  M.  Barnhart,  and 

Anthony  Ephremides* 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Washington,  DC  20375-5320 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NRL/MR/5521-95-7744 

9.  SPONSORINQ/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

OfEce  of  Naval  Research 

Arlington,  VA  22217 

10.  SPONSORINQ/MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

•University  of  Maryland  and  Kaman  Sciences  Coiporation 

12a.  DISTRIBUnON/AVAILABILfTY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  fAfejdmum  200  words) 

In  this  report,  we  consider  die  evaluation  of  data-packet  delay  in  wireless  integrated  voic^data  networks.  In  netwoiks  that  support 
circuit-switched  voice,  the  voice  occiq>ancy  process  satisfies  a  product-form  solution  under  reasonable  modeling  assumptions. 

Although  this  product-form  solution  provides  an  accurate  characterization  of  equilibrium  voice-trafHc  behavior,  it  does  not  directly 
provide  a  method  to  evaluate  data-packet  delay.  However,  examination  of  each  link  separately  in  a  manner  that  incorporates 
interaction  with  the  rest  of  the  network  permits  us  to  take  advantage  of  the  wireless  nature  of  the  networic  and  obtain  a  three-flow 
characterization  of  each  link,  which  also  satisfies  a  product-form  solution  and  is  hence  termed  a  ”mini-product-form"  solution.  By 
matching  the  values  of  these  flows  to  the  average  values  obtained  from  the  product-form  solution  of  the  entire  network,  we  obtain 
a  three-dimensional  Markov  chain  characterization  of  the  voice  occupancy  state  on  the  link,  which  permits  a  simpler  evaluation  of 
dau-packet  delay.  A  further  reduction  is  possible  by  converting  the  three-dimensional  chain  to  a  single-dimensional  one. 

Performance  results  demonstrate  that  these  models  provide  satisfactory  delay  estimates  that  also  appear  to  be  upper  bounds  on  delay. 

14.  SUBJECT  TERMS 

Conununications  network  Performance  evaluation 

Voice/data  integration  Product  form 

Radio  network 

15.  NUMBER  OF  PAGES 

37 

16.  PRICE  CODE 

17.  SECURfTY  CLASSFICATION 

OF  REPORT 

UNCLASSIFIED 

18.  SECURtTY  CLASSIFICATION 

OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

UNCLASSIFIED 

20.  LIMITATION  OF  ABSTRACT 

UL 

NSN  7540-01  -280-5500  Standard  Form  296  (Rav.  2-89) 

Praacriba  by  ANSI  Std  239-1 8 
298-102 


CONTENTS 


1.  INTRODUCTION . 1 

1 . 1  Earlier  studies  and  their  relation  to  this  report . 1 

1.2  Outline  of  the  Report . 2 

2.  INTEGRATED  WIRELESS  NETWORK  MODEL . 3 

2.1  Voice  traffic . 3 

2.2  Data  traffic . 4 

2.3  Data-packet  delay  evaluation  considerations  in  wireless  networks . 5 

3.  ON  THE  DEVELOPMENT  OF  DELAY  MODELS . 6 

3.1  An  exact  (but  impractical)  Markov  model . 6 

3.1.1  A  quasi-static  approximation . 9 

3.2  A  crude  data-traffic  model . 10 

4.  IMPROVED  APPROXIMATE  MODELS . 11 

4. 1  The  mini-product-form  model . 11 

4.1.1  Wired  vs  wireless  networks . 15 

4.2  A  reduced-load  model . 16 

5.  A  FURTHER  SIMPLMCATION  OF  THE  MODEL . 18 

6.  PERFORMANCE  RESULTS . 20 

6. 1  Accuracy  of  the  mini-product-form  model  characterization  of  the 

voice  process . 21 

6.2  Data-packet  delay  evaluation . 27 

6.2. 1  Data-packet  queue-size  distribution . 29 

7.  CONCLUSIONS . 31 

REFERENCES . . . 33 

Accesion  For  I 

NTIS  CRA&i  3 


OTIC  TAB  □ 

Unannounced  □ 

Justification 


By . . 

Distribution/ 


iii 

Availability  Codes 

Avail  and/or 

Dist 

Special 

i 

A-l 

Novel  Techniques  for  the  Analysis  of 
Wireless  Integrated  Voice/Data  Networks 


1.  Introduction 

The  accurate  prediction  of  average  data  delay  in  a  multihop  network  is  a  notoriously 
difficult  problem  that  has  been  extensively  addressed  in  the  literature  [1,2].  This  problem  is 
especially  difficult  for  integrated  voice/data  networks  in  which  the  quantity  of  resources 
available  for  data  traffic  (channels  in  the  case  of  wireline  networks  and  transceivers  in  the  case  of 
wireless  networks)  depends  on  the  current  state  of  the  voice  traffic.  Furthermore,  for  a  wireless 
network,  the  difficulty  is  compounded  by  the  interaction  amongst  the  transmissions  of 
neighboring  nodes. 

In  this  report  we  address  the  development  of  models  for  data-packet  delay  evaluation  in 
multihop  wireless  integrated  networks.  A  major  obstacle  to  the  development  of  accurate  delay 
models  has  been  the  lack  of  an  accurate  characterization  of  the  time-varying  nature  of  voice 
traffic  in  such  networks.  Such  a  characterization  is  necessary  because  models  based  simply  on 
the  expected  quantity  of  resources  available  for  data,  or  even  on  the  equilibrium  distribution  of 
the  quantity  of  resources  available  for  data,  can  be  highly  inaccurate.  Based  on  a  standard  model 
for  circuit-switched  voice  networks,  our  starting  point  is  the  use  of  the  product-form  solution  to 
characterize  the  state  of  voice-call  traffic.  Although  such  models  are  normally  used  to 
characterize  equilibrium,  rather  than  time-varying,  behavior,  we  show  that  the  product-form 
solution  can  be  used  in  a  novel  way  to  develop  an  approximate  Markovian  model  that  offers  a 
practical  and  useful  characterization  of  the  voice  occupancy  of  any  link  in  the  network,  which  in 
turn  provides  the  time-varying  “residual  capacity”  available  for  data,  and  hence  serves  as  the 
basis  for  the  evaluation  of  data-packet  delay.  We  compare  our  model  to  one  based  on  the  well- 
known  reduced-load  approximation,  and  find  significant  improvement  in  accuracy. 

1.1  Earlier  studies  and  their  relation  to  this  report 

This  report  is  part  of  our  continuing  study  of  voice/data  integration  in  wireless  networks, 
in  which  we  have  addressed  several  different  aspects  of  design  and  performance  analysis.  Some 
of  our  most  relevant  publications  are  summarized  below.  For  example,  in  [3-6]  we  developed 
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and  evaluated  the  WIMA  channel-access  protocol  for  integrated  networks.  The  Markovian 
system  model  developed  for  that  protocol  has  served  as  the  basis  for  the  multidimensional 
Markovian  model  for  system  evolution  in  integrated  wireless  networks  that  we  use  in  this  report 
for  data-packet  delay  evaluation.  In  [7-9]  we  addressed  the  problem  of  voice-call  admission 
control  in  integrated  wireless  networks,  taking  into  consideration  the  product-form  of  the 
equilibrium  voice-call  state  that  holds  under  reasonable  modeling  assumptions.  We  developed 
both  a  recursive  search  procedure  and  a  descent-search  technique  to  speed  up  the  determination 
of  the  optimal  admission-control  policy.  Our  use  of  the  product-form  solution  to  characterize 
equilibrium  behavior  in  those  studies  has  served  as  the  basis  for  its  use  in  dynamic  system 
models  in  the  current  report.  Also  of  major  significance  are  our  studies  of  ordinal  optimization 
methods  for  the  determination  of  good  admission-control  policies  in  integrated  wireless 
networks.  In  [10-13]  we  demonstrated  that  crude  analytical  or  simulation  models  with  low 
complexity  can  provide  nearly  optimal  performance,  although  the  delay  estimates  they  provide 
are  typically  poor.  In  the  present  report,  we  focus  on  the  development  of  models  for  data-packet 
delay  evaluation  in  integrated  wireless  networks.  Portions  of  this  work  have  been  presented 
recently  at  IEEE  INFOCOM’95  [14],  Many  of  our  recent  papers  are  available  by  anonymous  ftp 
from  mars.itd.nrl.navy.mil,  in  the  directory  “nrlpapers/code_5521/W.E.B.”,  or  on  the  WWW  at 
http://tang.itd.nrl.navy.mil/WEBindex.html. 

1.2  Outline  of  the  Report 

In  Section  2  we  discuss  our  integrated  voice/data  multihop  wireless  network  model,  and 
address  the  differences  between  wireless  and  wired  networks.  In  Section  3  we  discuss  issues  in 
the  development  of  delay  models  for  integrated  networks,  including  the  need  to  characterize  the 
time- varying  behavior  of  the  voice-call  process  (and  hence  the  residual  capacity)  over  individual 
network  links.  In  Section  4  we  develop  our  approximate  “mini-product-form”  model,  which 
provides  a  three-dimensional  first-order  Markov  chain  description  of  the  voice-call  process  in  the 
system;  this  provides  a  considerable  reduction  in  complexity  over  the  exact  Markov  chain 
description.  Then  in  Section  5  we  demonstrate  how  the  voice-call  description  can  be  reduced  to 
an  approximate  one-dimensional  first-order  Markov  chain,  with  little  loss  of  accuracy,  thus 
increasing  significantly  the  size  of  the  problems  that  can  be  handled  with  the  mini-product-form 
network  approach.  In  Section  6  we  present  performance  results  that  demonstrate  the  accuracy  of 
our  models,  and  in  Section  7  we  present  our  conclusions. 
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2.  Integrated  wireless  network  model 

In  wireless  networks  it  is  possible  to  establish  a  link  between  a  pair  of  neighboring  nodes, 
provided  that  each  has  a  transceiver  available  for  this  purpose.  Thus,  unlike  the  case  of  wired 
networks,  the  set  of  network  links  and  their  capacities  is  not  determined  a  priori,  but  can  adapt  to 
changing  traffic  demands.  This  is  a  crucial  difference  that  eventually  forms  the  basis  of  what  we 
propose  in  this  report.  We  assume  that  each  node  has  several  transceivers,  thereby  permitting 
flexibility  in  link  establishment,  while,  unfortunately,  also  complicating  performance  evaluation. 
In  this  report  we  do  not  address  the  protocol  issues  that  are  associated  with  call  setup  or  with  link 
activation,  but  instead  focus  on  the  development  of  a  mathematical  system  model  that  will  serve 
as  the  basis  for  data-packet  delay  evaluation  for  a  given  protocol  of  link  activation. 

We  consider  multihop  wireless  networks  in  which  circuit  switching  is  used  for  voice 
traffic  and  packet  switching  is  used  for  data  traffic.  After  describing  the  voice  and  data  traffic 
models,  we  address  in  more  detail  the  special  issues  associated  with  performance  evaluation  in 
wireless  networks. 

2.1  Voice  traffic 

We  consider  multihop  wireless  networks  in  which  voice-call  service  is  provided  by  the 
establishment  of  a  circuit  over  a  predetermined  path  between  the  originator  of  a  call  and  its 
destination  node.  This  requires  the  commitment  of  a  transceiver  at  every  node  along  the  path  for 
the  duration  of  the  call.  Let  T,-  denote  the  number  of  transceivers  at  node  i.  We  refer  to  calls  that 
use  circuit  j  <j  <  J)  as  calls  of  type  j;  they  arrive  to  circuit  j  according  to  a  Poisson  process 
with  rate  ?iy,  and  their  lengths  are  exponentially  distributed  with  parameter  |ij.  The  offered 
load  to  circuity  is  then  pj  =  /|xj.*  The  state  of  the  system  is  described  by  the  vector  jc  =  {xi, 

•••»  where  xj  is  the  number  of  calls  currently  active  on  circuit  j.  This  is  a  commonly  used 
circuit-switched  model  for  voice  networks  (see  [8, 9, 15, 16] ). 

A  central  controller  makes  the  decisions  on  whether  or  not  to  accept  calls  based  on 
perfect  knowledge  of  the  number  of  calls  of  each  type  that  are  currently  in  progress  (i.e.,  the 
system  state  jc),  and  hence  the  set  of  resources  that  are  available  for  new  calls.  The  transceivers 
needed  to  establish  a  circuit  are  acquired  simultaneously  when  the  call  arrives,  and  are  released 
simultaneously  when  the  call  is  completed.  Calls  are  blocked  when  one  or  more  nodes  along  the 
path  do  not  have  a  transceiver  available,  or  when  a  decision  is  made  not  to  accept  a  call  despite 
the  availability  of  transceivers. 

'  It  is  reasonable  to  have  values  of  pJ  greater  than  1  because  there  are  generally  several  transceivers  at  each  node. 
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If  we  restrict  the  voice  admission-control  policy  to  the  class  of  policies  known  as 
coordinate  convex  [17],  the  Markovian  nature  of  the  voice  traffic,  coupled  with  its  independence 
from  data  traffic,  results  in  a  product-form  solution  for  the  equilibrium  voice  state  of  the  system 
[18,  19].  Coordinate-convex  policies  are  a  reasonable,  easily  implemented  class  of  policies  [8, 
18].  Each  such  policy  is  specified  in  terms  of  the  set  of  admissible  states  Q.\  resources 
(transceivers)  freed  by  completed  calls  are  immediately  available  for  use  by  new  calls,  ongoing 
calls  do  not  get  rerouted,  and  a  new  call  is  admitted  if  and  only  if  the  state  to  be  entered  is  in  the 
admissible  region  (otherwise,  it  is  blocked  and  lost  from  the  system).  Although  the  optimal 
admission-control  policy  is  not  generally  coordinate  convex,  it  was  demonstrated  in  [18,  19]  for 
several  example  problems  that  the  best  coordinate-convex  policy  provides  nearly  optimal 
performance.  The  simplest  example  of  a  coordinate-convex  policy  is  the  example  of  an 
“uncontrolled”  system,  i.e.,  a  network  in  which  calls  are  admitted  as  long  as  resources  are 
available  to  support  them. 


The  equilibrium  state  probability  7i(-)  is  given  by 


7t(A:)  =  7c(0)n 


;=i 


AC  eQ 


where  7C(0)  =  (XareQ  11^=1  Ay  -)  normalization  constant  associated  with  the  admissible 

region  Q.  For  any  coordinate-convex  admission-control  policy,  it  is  straightforward  (though 
possibly  time  consuming)  to  evaluate  7C(0),  which  in  turn  permits  the  evaluation  of  performance 
measures  such  as  throughput  and  blocking  probability. 


2.2  Data  traffic 

Data  traffic  consists  of  single  fixed-length  packets  with  Poisson  arrival  statistics  (at  rate 
at  queue  k).  At  each  node  a  separate  queue  is  formed  for  packets  intended  for  each  of  its 
neighbors.  We  consider  only  the  single-hop  data  case,  i.e.,  relaying  is  not  addressed.^  As  in  the 
case  of  voice  traffic,  a  transceiver  is  needed  at  both  the  source  and  destination  nodes  throughout 
the  duration  of  a  data-packet  transmittal.  Voice  traffic  has  priority  over  data  traffic  in  the  sense 
that  the  decision  to  accept  a  voice  call  depends  only  on  the  number  of  voice  calls  (of  each  type) 


2  The  restriction  to  single-hop  data  traffic  is  not  really  a  major  limitation.  For  a  multihop  case,  we  would  derive  the 
data-packet  arrival  rate  at  each  queue  along  a  given  path  from  the  end-to-end  arrival  rates  (i.e.,  the  arrival  rates  of 
each  data-stream  type)  and  the  parameters  of  a  given  routing  algorithm.  Then  delay  over  multihop  data  paths,  for 
the  case  of  fixed  routing,  can  be  obtained  by  concatenating  the  delays  over  the  links  along  the  paths.  Of  course,  the 
implicit  assumption  of  independence  along  the  concatenated  links  makes  this  computation  an  approximate  one. 
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currently  in  progress  in  the  system,  and  is  independent  of  the  size  of  data  queues.  At  any  time 
instant,  the  transceivers  not  being  used  to  support  voice  traffic  are  available  to  support  data 
traffic.  Thus,  the  network  resources  available  for  data  depends  on  the  voice-call  statistics,  and  on 
the  voice-call  admission-control  policy. 

2.3  Data-packet  delay  evaluation  considerations  in  wireless  networks 

The  evaluation  of  data-packet  delay  is  complicated  by  the  flexible  manner  in  which  each 
node’s  transceivers  are  allocated  to  voice  and  data  traffic  in  wireless  networks.  Let  us  consider 
two  neighboring  nodes  a  and  b.  Typically,  at  an  arbitrary  time  instant,  some  of  each  node’s 
transceivers  will  be  supporting  voice  calls.  The  remainder  are  available  for  data  as  in  the 
movable-boundary  scheme  that  is  commonly  used  for  integrated  multiplexing  and  multiple 
access  [3, 20].  As  a  first  step  toward  developing  a  delay  model,  we  assume  that  any  transceivers 
not  occupied  by  voice  traffic  at  nodes  a  and  b  are  available  to  support  data  communication 
between  these  two  nodes.  We  thus  define  the  residual  capacity  Cab  available  for  data  on  “link 
(a,Z>)’’  as 

Cab  -  niin  {number  of  transceivers  not  occupied  by  voice  at  node  /}  . 

i=a,b 


In  general,  not  all  of  Cat  will  be  available  to  support  link  (a,b),  since  some  transceivers  at 
nodes  a  or  b  may  form  links  for  data  traffic  with  the  node’s  other  neighbors.  Typically,  a  link 
activation  [21-23],  frequency  management,  or  other  adaptive  demand-based  protocol  is  used  to 
partition  the  use  of  the  node’s  transceivers  among  all  the  links  that  the  node  may  form  with  its 
neighbors;  thus,  in  a  data-only  system,  our  link  (ajj)  would  receive  a  certain  fraction  (say  0)  of 
the  maximum  data  rate  that  can  be  supported  between  nodes  a  and  b.  In  the  integrated  network, 
one  possible  approach  would  be  for  link  (a,b)  to  receive  the  same  fraction  0  of  Cab-  To  avoid  the 
need  to  specify  a  particular  link-activation  schedule,  in  this  report  we  simply  assume  that  the 
entire  residual  capacity  of  the  nodes  (measured  in  terms  of  transceivers)  is  available  to  each  link 
of  interest.  It  would  be  straightforward  to  apply  our  approach  to  a  system  operating  under  any 
particular  link-activation  scheme^  by  using  this  reinterpretation  of  residual  capacity. 

In  Section  3.2  we  discuss  a  crude  model  for  delay  based  on  the  average  residual  capacity, 
and  in  Sections  4  and  5  we  develop  more-accurate  models  that  take  into  account  the  time  varying 
nature  of  the  residual  capacity. 


^  This  is  true  as  long  as  the  link-activation  scheme  does  not  depend  on  the  data-packet  queue  sizes  in  the  network. 


5 


3.  On  the  development  of  delay  models 

Although  performance  measures  such  as  voice  blocking  probability  and  throughput  can 
be  exactly  calculated  by  means  of  the  product-form  solution,  there  are  no  accurate  models 
available  for  calculating  data-packet  delay.  Accurate  models  require  a  characterization  of  the 
time-varying  behavior  of  the  voice  state.  Since  the  voice  state  of  the  system  is  Markovian,  it  is 
possible,  in  principle,  to  trace  the  distribution  of  the  voice  state  as  a  function  of  time,  and  thereby 
to  obtain  its  time-varying  behavior.  The  expected  data-packet  queue  size  (for  the  present  case  of 
single-hop  data  traffic  with  Poisson  arrivals)  can  then  be  studied  by  viewing  each  link  in  the 
system  as  a  variable  rate  queueing  system  with  service  rate  dependent  on  the  voice  process. 
Once  the  expected  data-packet  queue  size  is  obtained,  expected  delay  can  be  determined  by 
means  of  Little  ’  s  formula: 

,  ,  .  .  ,  ffdata  packet  queue  size} 

£{time  data  packet  spends  in  system}  =  - - — - 

where  is  the  data-packet  arrival  rate.  Unfortunately,  however,  such  a  Markov  model  for  voice 
traffic  would  be  practical  only  for  very  small  systems  (i.e.,  small  number  of  call  types  and  small 
maximum  number  of  calls  of  each  type).  In  addition,  the  study  of  the  behavior  of  a  queueing 
system  with  variable  rate  of  service  is  highly  nontrivial.  Thus  approximate  models  are  needed. 

In  Section  3.1  we  discuss  an  exact  Markov  model,  which  is  impractical  because  of  its 
dimensionality,  but  which  serves  as  the  basis  for  the  approximate  “mini-product-form”  model 
that  is  the  focus  of  this  report.  We  then  discuss  a  very  crude  M/D/1  approximation  in  Section 
3.2,  which  provides  remarkably  accurate  ordinal  rankings  of  control  policies,  but  poor  estimates 
of  actual  delay. 

3.1  An  exact  (but  impractical)  Markov  model 

In  this  section  we  outline  the  development  of  an  exact  Markov  model  for  the  numerical 
performance  evaluation  of  integrated  networks  of  the  type  described  in  Section  2.  We  make  use 
of  the  assumption  made  in  Section  2.3  that  the  entire  residual  capacity  of  any  pair  of  nodes  is 
available  to  the  link  that  connects  them.  Although  the  data-packet  delay  over  any  link  is 
independent  of  the  data  arrival  process  at  any  other  link  (because  we  are  considering  single-hop 
data  traffic),  the  delay  actually  depends  strongly  on  the  time-varying  residual  capacity  available 
for  data,  which  is  determined  by  the  voice  process  throughout  the  network.  Therefore,  a 
characterization  of  the  system  state  over  the  link  of  interest  must  incorporate  the  following 
quantities: 
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V(k)  =  {Vi(^),  V2(k),  Vj{k)},  where  Vj(k)  is  the  number  of  calls  of  type  j  that  are  in  the 

system  at  time  slot  k.  V(k)  determines  the  total  number  of  calls  of  all  types  that  are 
present  at  each  node  at  time  slot  k,  and  thus  determines  the  residual  capacity  of  any 
given  link  at  time  slot  k. 

R(k)  =  the  number  of  data  packets  in  queue  at  the  link  of  interest  at  time  slot  k.‘^ 

The  duration  of  a  time  slot  is  the  time  required  for  a  transceiver  to  transmit  one  fixed-length  data 
packet  (all  transceivers  transmit  at  the  same  rate). 

The  analysis  of  our  integrated  network  is  similar  to  that  developed  for  the  movable¬ 
boundary  Wireless  Integrated  Multiple  Access  (WIMA)  protocol  [3],  and  can  proceed,  in 
principle,  as  follows.  The  behavior  of  the  link  of  interest  can  be  characterized  by  the  Markov 
chain  {R(k),  V(^)},  which  has  transition  probabilities  Pr{/?(A:+1),  V(k+1)  I  R{k),  V(k)}.  The 
steady-state  distribution  of  the  system  state  can  then  be  obtained  by  repeated  iteration  of  the 
transition  probability  matrix  until  equilibrium  is  reached.  In  the  numerical  evaluation  it  is 
necessary  to  truncate  the  data-packet  queue  size  (which  is  inherently  unbounded)  at  a  finite 
value,  which  should  be  chosen  sufficiently  large  so  that  truncation  effects  on  the  resulting 
equilibrium  probability  distribution  are  negligible;  no  such  truncation  is  needed  for  the  voice-call 
process,  which  is  finite  because  it  is  limited  by  the  number  of  transceivers  at  the  nodes  or 
possibly  by  an  admission-control  policy.  The  expected  data-packet  delay  at  this  link  can  be 
evaluated  by  determining  the  expected  data-packet  queue  size  in  equilibrium.  The  numerical 
evaluation  outlined  here  is  repeated  at  each  link  for  which  expected  delay  values  are  desired.  It 
is  a  complex  and  cumbersome  evaluation,  and  is  presented  here  just  to  sketch  out  the  process  by 
which  an  exact  evaluation  could  be  carried  out  if  the  size  of  the  network  were  small. 

To  reduce  the  dimensionality  of  the  transition-probability  matrix  (which  may  be 
prohibitively  large),  we  observe  as  in  [3]  that  the  data  traffic  depends  on  the  voice  state,  whereas 
voice  is  independent  of  data.  Therefore,  we  have 

PT{R(k  +  l),V(k  +  l)\R(k),Vik)}  =  Pr {R(k ■¥  1) \  R(k),V(k)}  Pr{y(/fc-H)l  V(k))  . 

The  transition  from  R(k)  to  Rik+l)  depends  on  V{k)  (because  V(k)  determines  the  residual 
capacity),  but  not  on  V{k+l).  The  transition  from  V(k)  to  V(k+l)  does  not  depend  on  R(k)  or 
R{k+1). 


*  A  single  conceptual  queue  at  link  (a,b)  is  used  to  represent  all  data  traffic  in  either  direction  between  nodes  a  and 
b.  This  entails  no  loss  of  generality  since  we  assume  that  a  full  transceiver  is  required  at  both  the  transmitting  and 
receiving  nodes  of  a  data  link. 
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We  observe  that  the  transitions  from  slot  k  to  slot  ^+1  can  be  modeled  as  a  two-step 
process.  Data  transitions  are  considered  first.  Given  Fr{R(k),  V{k)},  we  first  determine 
Pr{R(k+l),  V'(^)}.  The  evolution  of  the  data-packet  process  between  slots  k  and  ^-i-l  can  be 
described  as  follows: 


R(k  +  l)  =  [/?(A:)-c(it)f +  A(jt), 

where  c{k)  is  the  residual  capacity  of  the  link  of  interest  in  slot  k  (determined  directly  from  F(/:)), 
A(k)  is  the  number  of  data-packet  arrivals  in  slot  k  (these  packets  cannot  be  transmitted  before 
slot  it+1),  and  [jr]+  =  max{j:,  0}.  We  assume  that  data-packet  arrivals  form  a  Bernoulli  process 
with  at  most  one  packet  per  slot^  with  an  arrival  probability  of  thus 

Pr{yl(it)  =  l)  =  Xj  and  Pr{A(it)  =  0}  =  l-Xj. 

Therefore,  the  data-packet  transition  probabilities  (conditioned  on  the  residual  capacity)  can  be 
directly  written  as: 


Pr{/?(it-l-l)  =  I/?(jt)-c(/:)ric(it)}  =  l-Xd, 


and 


Pr[R(k  +  \)  =  [R{k)-cik)Y  +  l\c(k)]  =  X^, 

from  which  we  obtain  Pr{/?(A:-t-l)  =  nl/?(fc)  =  i,V(^)  =  v}.  The  joint  distribution  Pr{/?(^-t-l), 
V(k)}  is  then  determined  as  follows: 

n+T 

Pr{/?(^-H)  =  n,V(it)  =  v}  =  '^?T{R(k  +  l)  =  n\R(k)  =  i,V(k)  =  v}  PT{R(k)  ==  i,Vik)  =  v) , 

i=n-l 


where  T  is  the  maximum  residual-capacity  value  that  the  link  of  interest  can  have  (namely  the 
maximum  of  the  number  of  transceivers  at  the  two  nodes  forming  this  link). 

Next,  the  voice  transitions  are  considered.  Given  Pr{/?(fc-«-l),  V(k)}  we  determine 
PT{R(k+l),  V(ik+1)}.  Since  the  voice  transitions  are  independent  of  the  data  traffic,  the  same 
transition  probability  matrix  is  used  for  all  values  of  R(k+l).  Thus,  the  following  is  evaluated: 


5  Other  distributions  for  could  certainly  be  considered  as  well,  e.g.,  a  Poisson  arrival  distribution  (truncated  at 
some  finite  value  to  facilitate  computation). 
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V 

^  max 

^x{R{k  + 1)  =  j,V{k  + 1)  =  H'}  =  ^Pr{V(k  +  l)  =  w\  V(k)  =  v}  FT{R(k  + 1)  =  y,y(jt)  =  v}, 

v=0 

where  the  summation  is  actually  a  7-dimensional  summation  over  all  possible  voice  states,  which 
range  from  the  empty  voice  state  (v  =  0)  to  voice  states  at  the  boundary  of  the  admissible  region 
(denoted  loosely  as  Vmax)-  The  equilibrium  distribution  of  the  system  state  is  determined  by 
repeating  this  two-step  iteration  until  convergence  is  achieved,  i.e.,  until  there  is  negligible 
change  in  the  distribution,^ 

Although  this  model  is  exact  (except  for  the  assumption  that  the  entire  residual  capacity 
is  available  for  the  link  of  interest),  its  use  is  impractical  because  of  the  high  dimensionality, 
which  results  from  the  need  to  model  the  7-dimensional  voice  state. 

3.1.1  A  quasi-static  approximation 

To  reduce  the  complexity  of  this  evaluation,  it  is  possible  to  consider  an  approximate 
quasi-static  model  that  exploits  the  fact  that  the  voice-call  process  changes  much  more  slowly 
than  the  data-packet  process;  e.g.,  in  practical  systems  the  voice-call  duration  may  typically  be 
several  orders  of  magnitude  greater  than  the  data-packet  time  slot.  First,  we  assume  that  V(k)  is 
constant  for  all  time.  For  each  possible  value  of  V(k)  we  can  then  determine  the  equilibrium 
distribution  of  R(k)  by  numerical  iteration  of  the  transition  probability  matrix  corresponding  to 
the  value  of  residual  capacity  associated  with  V(k).  From  this  equilibrium  distribution  we  can 
determine  the  expected  queue  size,  and  hence  expected  delay,  associated  with  each  voice-call 
state  V(k).  Finally,  we  can  average  the  distribution  of  R(k)  over  all  possible  values  of  Vik), 
whose  distribution  can  be  determined  exactly  by  means  of  the  product-form  distribution  without 
requiring  the  iterative  approach  discussed  here.  However,  unless  the  instantaneous  residual 
capacity  in  all  voice  states  is  greater  than  the  offered  data  load,  the  resulting  expected  delay  will 
be  infinite;  e.g.,  if  any  state  with  nonzero  probability  provides  a  residual  capacity  of  zero  on  the 
link  of  interest,  any  finite  data  rate  on  that  link  will  produce  infinite  expected  delay.'^  Our  mini¬ 
product-form  model,  which  incorporates  dynamic  system  behavior,  only  requires  that  the  offered 
data  traffic  be  less  than  the  expected  residual  capacity,  and  does  not  require  that  any  transceivers 
be  reserved  for  the  exclusive  use  of  data. 


®  In  our  studies  of  the  mini-product-form  model,  we  have  stopped  the  iteration  process  when  the  change  in  expected 
queue  size  between  iterations  n  and  n  +  10  was  less  than  10'^  percent. 

’  This  approach  was  used  to  evaluate  data-packet  delay  for  the  WIMA  protocol  in  [5]  for  “movable-boundary” 
examples  in  which  some  resources  were  reserved  for  exclusive  use  by  data  traffic. 
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3.2  A  crude  data-traffic  model 


For  each  voice-call  admission-control  policy,  we  can  compute  the  expected  residual 
capacity  at  every  link  by  means  of  the  product-form  solution.  If  we  make  the  obviously 
inaccurate  approximation  that  the  data  packets  at  link  {a,b)  are  served  at  the  constant  rate  of  the 
average  residual  capacity  of  the  link,  Cat,  the  resulting  system  is  an  M/D/1  queueing  system  in 
which  the  delay  (not  including  service  time)  is 


2)1(1 -p)’ 

Here,  p  =  X/|i,  where  X  is  the  Poisson  data-packet  arrival  rate  and  p  =  l/(time-slot  duration)  = 
Cab- 


This  delay  model  is  deficient  in  several  ways.  Most  significant  is  the  implicit  assumption 
that  p,  the  data-packet  service  rate  at  link  {a,b),  is  constant  at  the  expected  residual-capacity 
value  Cab  (the  fact  that  Cab  is  not,  in  general,  integer-valued  is  a  minor  point  here),  whereas  the 
number  of  transceivers  available  in  the  real  system  varies,  based  on  the  voice  state.  Thus  data- 
packet  queue  sizes  may  increase  rapidly  when  the  voice-call  occupancy  is  high  and  decrease 
rapidly  when  it  is  low.  This  model  also  ignores  the  duration  of  voice  calls®  and  hence  the  statis¬ 
tics  of  the  amount  of  time  spent  in  each  voice-occupancy  system  state.  Despite  the  poor  predic¬ 
tive  quality  of  this  delay  estimate  (see  Section  6.2),  we  have  found  that  this  model  performs 
surprisingly  well  as  a  method  for  ordinal  optimization,  i.e.,  for  the  ranking  of  the  voice-call 
admission-control  policies  with  respect  to  the  achieved  average  data  delay  computed  under  this 
model;  in  fact,  the  resultant  ranking  was  in  almost  perfect  agreement  with  the  ranking  found  by 
extensive  accurate  simulations  (which  do  not  make  the  M/D/1  approximation)  [12].  However, 
our  objective  in  the  present  report  is  to  improve  the  quality  of  the  actual  delay  estimate. 

The  high  complexity  of  the  exact  Markovian  model  and  the  inaccuracy  of  the 
approximate  models  discussed  in  this  section  point  to  the  need  for  the  development  of  accurate 
and  practical  models  for  delay  performance.  Such  models  must  address  the  time-varying  nature 
of  the  residual  capacity,  without  being  too  complex  for  feasible  computations.  In  Section  4  we 
discuss  the  development  of  such  a  model. 


®  It  uses  only  the  ratio  of  the  voice  arrival  rate  to  its  service  rate,  and  is  independent  of  the  actual  arrival  rate  and 
service  rate  values. 
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4.  Improved  approximate  models 

The  use  of  the  exact  Markovian  model  discussed  in  Section  3.1  is  impractical  because  of 
its  high  dimensionality,  which  results  from  the  need  to  track  the  voice-call  state  throughout  the 
network  while  evaluating  data-packet  delay  at  a  single  link  of  interest.  The  system  state  at  the 
link  of  interest  is  given  by  V(^)},  in  which  the  voice  state  V(^)  (for  a  system  with  J  voice 

circuits)  is  a  7-dimensional  first-order  Markov  chain.  In  this  section  we  discuss  an  approximate 
model  in  which  a  three-dimensional  first-order  Markov  chain,  which  we  refer  to  as  the  “mini- 
product-form”  model,  is  used  to  describe  the  voice-call  state  at  any  link  in  the  network,  resulting 
in  a  practical  methodology  for  performance  evaluation.  Our  model  is  similar  in  principle  to  the 
well-known  reduced-load  approximate  models  [15,  16,  24-28],  which  we  discuss  in  Section  4.2. 
However,  we  demonstrate  in  Section  6  that  our  mini-product-form  model  is  generally  more 
accurate  than  a  time-varying  model  based  on  the  reduced-load  approximation. 

4.1  The  mini'product-form  model 

Our  objective  is  to  find  a  description  of  the  voice-call  process  on  the  link  of  interest  that 
is  simpler  than  that  provided  by  the  7-dimensional  V{k),  but  yet  maintains  a  high  degree  of 
accuracy.  In  doing  so,  we  have  developed  an  approximate  model  in  which  we  consider  each  link 
separately,  in  terms  of  both  voice  and  data  traffic.  The  objective  is  to  characterize  the  time- 
varying  capacity  available  for  data  traffic  over  each  such  link,  while  taking  into  account  the 
effects  of  voice  traffic  throughout  the  entire  network.  In  particular,  the  voice  traffic  carried  by 
the  two  nodes  that  form  the  link  of  interest  (nodes  a  and  b)  must  be  modeled.  The  traffic 
supported  by  these  nodes  is  the  superposition  of  the  traffic  carried  by  all  circuits  passing  through 
them.  The  throughput  actually  delivered  by  each  of  these  circuits  is  generally  less  than  the 
offered  traffic  level  because  of  blocking,  both  at  the  nodes  of  interest  and  at  other  nodes.  In 
recent  years  a  number  of  “reduced-load”  approximations  have  been  developed  to  account  for 
blocking  at  other  links  or  nodes  in  the  network.  These  are  based  on  the  assumption  that  blocking 
occurs  independently  at  each  node,  an  assumption  that  is  asymptotically  accurate  for  large 
networks  with  a  large  number  of  transceivers  per  node  and  a  large  number  of  call  types,  where 
the  ratio  of  traffic  to  capacity  is  kept  fixed.  We  discuss  the  application  of  reduced-load 
approximations  to  our  problem  in  Section  4.2.  But  first  we  develop  our  mini-product-form 
model,  which  is  the  main  contribution  of  this  report. 

Consider  a  link  defined  by  node  pair  ab,  as  shown  in  Fig.  1,  where  nodes  a  and  b  have  Ta 
and  Tb  transceivers,  respectively.  The  flexibility  of  wireless  communications,  in  the  sense  that  a 
transceiver  can  form  a  link  with  any  of  its  neighbors  that  also  have  an  available  transceiver. 
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eliminates  the  need  to  assign  a  priori  any  particular  neighbor  or  call  type  to  a  transceiver.  We 
define  three  voice-traffic  flows  that  involve  link  {a,b)  as  follows: 

fab  =  flow  passing  through  both  nodes  a  and  b  (i.e.,  the  sum  of  the  expected  number  of 
calls  on  all  circuits  traversing  link  {a,b),  whether  terminating  or  starting  at  either 
node  a  or  Z?  or  whether  continuing  through  either  node  in  the  in-bound  or  out-bound 
direction), 

fa  =  flow  passing  through  node  a  but  not  node  b  (i.e.,  the  sum  of  the  expected  number  of 
calls  on  all  circuits  traversing  node  a,  whether  terminating,  starting,  or  continuing 
through  that  node  but  not  involving  node  b), 

fb  =  flow  passing  through  node  b  but  not  node  a. 


Fig.  1  —  The  “mini-product-form”  network. 


These  flows  would  be  artificial  and  useless  in  a  wired  network.  They  are  very  real  and 
useful,  however,  in  the  context  of  a  wireless  network  because  they  provide  a  way  to  distinguish 
between  the  time-average  number  of  transceivers  at  each  node  dedicated  to  service  link  (a,b)  and 
those  assigned  to  service  calls  that  simply  go  through  nodes  a  orb  without  traversing  link  {a,b). 
Note  that  the  total  flow  through  node  a  is  then/^i,  +fa,  and  the  expected  residual  capacity  at  node 
a  is  Ta  —  (fab  "^fa)-  Thus  we  can  consider  a  three-dimensional  voice  state  description  of  these  two 
nodes  at  any  time  instant  given  by  a:  =  [xa,  Xb,  Xab},  where 

Xab  =  number  of  voice  calls  currently  passing  through  both  nodes  a  and  b, 

Xa  =  number  of  voice  calls  currently  passing  through  node  a  but  not  node  b, 

Xb  =  number  of  voice  calls  currently  passing  through  node  b  but  not  node  a. 

Note  that  the  total  number  of  calls  currently  passing  through  node  a  is  +  Xa,  and  the 
instantaneous  residual  capacity  of  node  a  isTa-  (Xab  +  Xa). 

Unfortunately,  the  random  process  x  is  not  Markovian.  The  arrival  processes  associated 
with  the  three  flows  are  not  Poisson  because  of  blocking  in  the  network.  However,  we  make  the 
approximation  that  these  flows  are  characterized  by  Poisson  arrival  processes,  at  rates  chosen  to 
match  the  traffic  levels  in  the  overall  network.  In  doing  so  we  are,  in  effect,  approximating  the 
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distribution  of  the  interarrival  times  for  the  voice  traffic  on  a  given  link  by  its  maximum  entropy 
estimate  [29].  As  in  the  overall  network,  calls  may  be  blocked  because  of  a  lack  of  resources  or 
because  of  the  control  policy.  This  two-node  three-flow  system  (under  the  Poisson  arrival 
assumption)  is  also  a  product-form  system,  and  we  refer  to  it  as  our  “mini-product-form”  model. 
By  focusing  on  the  immediate  neighborhood  of  two  nodes,  we  attempt  to  capture  some  of  the 
interactions  that  complicate  the  network  behavior. 

We  conjecture  that  the  behavior  of  the  mini-product-form  network  associated  with  nodes 
a  and  b  closely  approximates  the  behavior  of  nodes  a  and  b  in  the  overall  network,  provided  that 
appropriate  system  parameters  are  chosen.  In  particular,  our  objective  is  to  find  the  offered 
traffic  levels  p' =  {Pa,Pi„Pai,}  for  the  mini-product-form  network  that  produce  the  same 
expected  traffic  flows  (  /  =  {faifb'fab^)  those  of  the  product-form  solution  for  the  network  as 
a  whole  operating  under  the  offered  load  p  =  {pi,  ...,p7}.  Here  we  discuss  the  simple  iterative 
procedure  we  have  developed  to  determine  the  offered  load  p',  and  in  Section  6  we  demonstrate 
that  the  mini-product-form  does,  indeed,  provide  a  good  characterization  of  system  behavior. 

Flow-Matching  Algorithm: 

Guess  initial  values  for  p'. 

Repeat  until  (error  <  e): 

{Calculate  via  product  form  the  flows  f^,  f{,,  and  fab  in  the  mini-product-form 
network  with  offered  load  p'; 
error  =  (/^  -  /^'  -H  /^  -  -  f'^)/3; 

adjust  p'  proportional  to  individual  errors:  fa- fa,  fb-fb,  and  fab- fab,  ie-. 
Pa  •=Pa +5{/a  ~/a)’  Pb  ■=Pb'^^{fb  ~  fb),  Pab '-Pab'^^ifab  ~  fab)-  ) 

The  primed  flows  f^,  fb,  and  f^  denote  the  flows  in  the  mini-product-form  network 
operating  under  the  offered  load  p^,  p^,  and  p^b-  We  have  found  that  reasonable  initial  values 
for  Pa,  Pb,  and  Pab  are  the  desired  flows /a, /t,  and/af,.  In  our  experiments,  we  have  used  a 
proportionality  coefficient  of  5  =1.0  in  adjusting  p',  and  have  typically  used  an  error  limit  value 
ofe=  10-6.9 

This  algorithm  converges  rapidly  to  provide  accurate  values  for  p'.  For  example.  Figure 
2  shows  the  number  of  iterations  that  were  required  for  our  iterative  procedure  to  converge;  the 
example  considered  is  that  of  link  (5,13)  in  the  network  of  Fig.  3,  which  will  be  discussed  in 
Section  6.  The  figure  shows  that  the  number  of  iterations  required  to  converge  to  within  e  of  the 


’Use  of  e  =  10"^  results  in  values  of  p'  that  differ  from  those  for  e  =  10^  only  in  the  third  or  fourth  decimal  place. 
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desired  flows  increases  with  offered  load,  but  is  generally  not  prohibitive.  Although  the  number 
of  iterations  depends  on  the  offered  load,  the  time  required  to  perform  a  single  iteration  of  this 
procedure,  i.e.,  to  find  the  product-form  solution  for  a  single  parameter  set  and  adjust  p',  depends 
primarily  on  the  number  of  transceivers,  which  in  turn  determine  the  number  of  states  that  must 
be  evaluated.  For  example,  with  four  transceivers  per  node  there  are  55  states,  and  a  single 
iteration  of  this  procedure  take  about  5.5  ms;  with  eight  transceivers  per  node  there  are  285 
states,  and  an  iteration  takes  about  14.4  ms;  and  with  16  transceivers  per  node  each  iteration 
requires  nearly  81  ms  to  evaluate  the  1785  states  associated  with  the  mini-product-form  model. 
The  rate  of  growth  in  the  time  required  to  perform  an  iteration  reflects  the  exponential  growth  of 
the  state  space  with  increasing  numbers  of  transceivers  per  node.  In  Section  5,  we  show  that  the 
number  of  states  in  the  mini-product-form  model  is  proportional  to  the  cube  of  the  number  of 
transceivers  at  a  node  (assuming  both  nodes  have  an  equal  number  of  transceivers).  Hence,  this 
approach  is  feasible  for  the  moderately  sized  examples  used  in  this  report.  However,  for  larger 
problems  an  approximation,  such  as  the  reduced-load  approximation  (see  Section  4.2),  will  have 
to  be  used. 


Fig.  2  —  Iterative  procedure  convergence  vs  offered  load  p  . 

Thus  although  the  product-form  solution  is  normally  used  to  characterize  equilibrium 
behavior,  we  use  it  as  the  basis  for  our  dynamic  time-varying  model.  By  matching  the  traffic 
flows,  we  attempt  to  capture  some  of  the  effects  of  blocking  at  remote  nodes.  Once  the  values  of 
p' ,  p^,  and  pQ^  are  determined,  the  time- varying  behavior  of  the  voice-call  link  occupancy  can 
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be  determined  by  numerical  solution  of  the  resulting  three-dimensional  Markov  chain  for  the 
mini-product- form  network  by  following  the  procedure  outlined  in  Section  3.1  for  the  overall 
network.  Doing  so  is  considerably  easier  than  solving  the  J-dimensional  Markov  chain  for  the 
overall  network,  particularly  when  the  number  of  call  types  is  large. 

A  potential  limitation  of  the  mini-product-form  model  is  that  its  use  requires  the 
evaluation  of  the  normalization  constant  associated  with  the  product-form  solution,  which  is 
computationally  intensive  for  large  networks.  However,  it  is  anticipated  that  an  acceptable 
estimate  of  the  voice-call  state,  from  which  the  parameters  for  our  mini-product-form  model  can 
be  obtained,  can  be  evaluated  by  using  a  reduced-load  approximation  to  obtain  the  equilibrium 
state  description.  In  Section  4.2  we  discuss  the  application  of  the  reduced-load  approximation  in 
a  somewhat  unconventional  manner  to  obtain  an  alternative  approximate  time-varying  model  for 
our  network,  but  show  that  its  performance  is  inferior  to  that  of  our  model.  In  Section  6  we 
discuss  performance  results  based  on  the  mini-product-form  model,  and  demonstrate  its  accuracy 
over  a  wide  range  of  situations. 


4.1.1  Wired  vs  wireless  networks 


In  developing  the  mini-product-form  model  for  wireless  networks  we  have  been  able  to 
exploit  the  flexibility  of  wireless  communication,  which  permits  a  transceiver  to  form  a  link  with 
any  of  its  neighboring  nodes.  This  flexibility  has  permitted  the  characterization  of  the  voice 
traffic  over  the  link  of  interest  in  terms  of  a  three-flow  process.  In  contrast,  the  development  of 
accurate  models  for  wired  networks  is  more  difficult  because  of  the  fixed  nature  of  the  network 
link  capacities,  as  we  now  explain. 


In  an  attempt  to  obtain  an  accurate  model  for  voice-call  traffic  over  individual  links  in 
product-form  wired  networks,  we  investigated  the  use  of  an  M/MJm/m  queueing  model,*®  which 
is  commonly  used  to  model  voice  systems.  Here,  the  parameter  m  is  simply  the  number  of 
channels  supported  by  the  link  of  interest.  The  throughput  of  an  M/M/m/m  queue  is 


ThroughputMAi//n/m  =  P 


1- 


p^/m! 


^n=0^ 


In  a  manner  analogous  to  that  used  in  our  wireless  network  model,  we  modeled  an  arbitrary  link 


An  M/M/m/m  queue  is  one  in  which  both  the  airival  and  service  times  are  exponential,  the  number  of  servers  is  m, 
and  the  total  number  of  customers  permitted  in  the  system  is  also  m.  Thus  this  is  a  loss  system  in  which  queueing  is 
not  permitted;  customers  are  rejected  if  they  arrive  when  all  servers  are  busy. 


15 


as  an  M/M/m/m  queue,  and  we  determine  the  offered  load  p'  (a  single  parameter  in  this  case)  so 
as  to  reflect  the  interactions  of  the  voice  traffic  on  the  link  of  interest  with  the  rest  of  the  network. 
Specifically,  p'  was  chosen  such  that  the  expected  voice-call  throughput  supported  by  the  link 
(in  isolation,  using  the  M/M/m/m  model)  was  equal  to  that  which  is  generated  by  the  overall 
product-form  solution,  which  characterizes  global  network  performance.  Although  this  model 
performed  well  in  some  cases,  errors  were  typically  much  larger  than  those  obtained  by  using  the 
mini-product-form  solution  for  wireless  networks. 

It  is  not  surprising  that  the  mini-product-form  solution  for  wireless  networks  performs 
better  than  the  M/M/m/m  model  for  wired  networks,  because  the  mini-product-form  solution 
matches  not  only  the  flow  across  the  link  of  interest,  but  also  the  adjacent  flows.  To  match  flows 
in  a  wired  network  would  require  a  separate  flow  for  each  link  adjacent  to  the  link  of  interest.  In 
some  cases  this  might  require  more  than  J  (the  number  of  voice  circuits)  flows,  thus  resulting  in  a 
model  with  complexity  greater  than  that  of  the  Markovian  model  for  the  overall  network  (which 
has  one  flow  per  circuit  type).  The  feature  of  the  wireless  network  that  permits  a  reasonably 
accurate  system  description  using  three  flows  is  the  fact  that  the  transceivers  at  a  node  are  not  a 
priori  assigned  to  specific  links,  but  may  form  links  with  any  of  their  neighboring  nodes.  Hence, 
a  single  flow  can  incorporate  the  effects  of  all  flows  adjacent  to  the  link  of  interest. 

4.2  A  reduced-load  model 

The  “reduced-load”  model  has  been  proposed  for  the  evaluation  of  traffic  statistics  at 
individual  links  or  nodes  in  multihop  circuit-switched  product-form  networks  [15,  16,  24-28].  If 
there  were  no  blocking  in  the  network,  the  average  load  offered  to  a  node  would  be  simply  the 
sum  of  the  offered  load  on  each  circuit  that  uses  the  node.  That  is,  Pj  >  where  Pa  = 

{circuits  that  use  node  a],  pj  =  is  the  offered  load  on  circuit  j  and  is  the  average  load 
offered  to  node  a.  However,  as  a  result  of  resource  limitations  (e.g.,  a  finite  number  of 
transceivers)  or  admission-control  policies,  the  load  offered  to  a  node  can  be  “thinned”  or 
reduced  by  blocking  at  other  nodes.  The  reduced-load  approximation  estimates  the  thinned 
average  offered  load  at  a  node  in  such  a  network  by  making  the  assumption  that  blocking  occurs 
independently  at  each  node  so  that  the  estimated  offered  load  at  node  a  is 

Pi  =  Zpj 

yePfl  /ey-a 

where  the  product  is  taken  over  all  nodes  I  that  are  elements  of  circuit  j  except  node  a,  and  P^(l) 
is  the  probability  of  blocking  at  node  /.  Thus,  the  product  term  represents  the  probability  of  call 
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acceptance.  The  remaining  difficulty  is  calculating  (/).  Typically  a  method  of  repeated 
substitution  is  used  to  solve  the  Erlang  loss  formula 


PbH) 


where  Ti  is  the  number  of  transceivers  at  node  1. 


For  comparison  purposes,  we  have  used  the  reduced-load  approximation  to  estimate  the 
load  that  should  be  offered  to  a  node  in  isolation  to  yield  the  same  statistics  as  would  have 
been  seen  at  node  a  operating  as  an  element  of  the  whole  network.  For  relatively  small 
networks,  we  can  calculate  the  product-form  solution  for  the  network  to  obtain  exact  values  of 
circuit  blocking  probability  (denoted  If  (j)  for  circuit  j)  and  node  blocking  probability  (denoted 
Pb(o)  for  node  a),  thus  eliminating  the  need  to  use  repeated  substitution."  Based  on  the 
independence  assumption  used  in  this  method,  we  observe  that 


1-PtO) 

!-/?(<!)■ 


We  then  use  Eq.  (1)  to  obtain 


Pa 


1 


Following  that,  we  can  compute  the  state  occupancy  distribution  Ti(ya)  for  node  a  in  isolation 
since  it  has  the  product-form  given  by 


(Pa)”/ 
^n=0  /n\ 


=o,i,...,r^. 


where  y^,  the  total  number  of  calls  currently  passing  through  node  a,  is  equal  to  Xa  +  Xat  in  the 
mini-product-form  model,  as  defined  earlier.  *2 


' '  For  examples  too  large  to  compute  the  normalization  constant  (in  which  case  the  reduced-load  approximation 
generally  performs  well),  the  standard  method  of  repeated  substitution  would  be  necessary.  In  fact,  this  method 
could  be  used  to  obtain  the  parameters  for  our  mini-product-form  model  in  large  examples. 

12  Our  use  of  the  reduced-load  model  is  unconventional  in  that  we  assume  that  the  normalization  constant  is  known. 
Usually  the  reduced-load  model  is  used  to  approximate  the  normalization  constant,  as  discussed  in  Section  4.1. 
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To  obtain  our  reduced-load  estimate  of  the  statistics  of  link  {a,b)  in  a  wireless  network, 
we  make  the  further  approximation  that  the  state  occupancy  distributions  niy'a)  and  K{yb)  are 
independent,  i.e.,  that  K{ya,  yt)  =  Tt(ya)  The  residual  capacity  available  for  data  on  link 

(a,b)  is  Cah  =  min{(ra  -  ja),  (Tb-yt)},  and  its  distribution  can  be  obtained  from  the  approximate 
joint  node  occupancy  distribution  n{ya,  yt)- 

In  summary,  in  our  reduced-load  approximation  the  voice  traffic  over  link  (a,b)  is 
modeled  as  the  superposition  of  two  voice  flows  at  nodes  a  and  b,  which  are  conditionally 
independent,  given  the  offered  load  to  the  network.  By  contrast,  in  the  mini-product-form  model 
the  correlation  in  the  flows  through  nodes  a  and  b  is  reflected  by  the  flow  fab-  Taking  this 
correlation  into  consideration  results  in  a  much  better  estimate  of  system  performance,  as  shown 
in  Section  6. 

5.  A  further  simplification  of  the  model 

Although  the  three-dimensional  mini-product-form  model  for  voice  traffic  results  in  a 
substantial  reduction  in  model  complexity  as  compared  with  the  exact  7-dimensional  model,  the 
complexity  remains  high.  In  this  section  we  discuss  the  use  of  a  one-dimensional  Markov  model 
for  voice  traffic  over  a  link,  whose  use  introduces  only  small  errors  as  compared  to  the  three- 
dimensional  mini-product-form  model. 

Let  us  again  consider  the  mini-product-form  network  shown  in  Fig.  1,  and  in  particular 
the  case  in  which  Ta  =  Tb  =  T  (i.e.,  the  number  of  transceivers  at  both  nodes  a  and  b  is  T).  The 
number  of  voice  states  can  be  evaluated  as  follows.  For  a  given  value  of  Xab,  the  number  of 
possible  values  of  jc^  is  (T  -  Xab  +1)  because  the  number  of  transceivers  occupied  by  voice  at 
node  a  can  range  from  0  to  J;  the  same  is  true  for  the  number  of  possible  values  of  Xb-  Thus  the 
number  of  possible  [Xa,  Xb)  pairs  for  the  specified  value  of  Xab  is  (T  -  Xab  +  1)^-  Summing  over 
all  possible  values  of  Xab  to  find  the  total  number  of  states  of  the  form  [Xa,  Xb,  Xab)  yields; 

T’+l 

number  of  states  in  mini-product-form  model  = 

n=l 

^  (T+l)(T+2)(2T+3)  ^ 

Overall,  the  total  number  of  states  possible  for  this  link  is  this  number  multiplied  by  D^ax  +  1. 

This  independence  approximation  can  result  in  errors  in  the  joint  distribution  that  are  significantly  larger  than 
those  based  on  the  reduced-load  approximation  at  the  individual  nodes. 
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where  Dmax  is  the  value  at  which  the  data-packet  queue  size  is  truncated  in  the  numerical 
evaluation.  At  high  throughput  levels  we  have  used  values  of  Dmax  as  great  as  1000  to  maintain 
the  accuracy  of  the  calculations.  This  dimensionality  limits  the  applicability  of  the  mini-product- 
form  solution  to  relatively  small  problems. 

We  make  the  observation  that  the  data-packet  queueing  process  on  link  {a,b)  depends  on 
the  voice-call  process  only  through  the  residual  capacity  Cab  available  for  data,  which  at  any 
given  time  is 

Cab  =  nun{[ Ta-iXa+Xab)]  ,[Tb-(Xb+Xab)]},  (2) 

where  the  time  dependence  of  Cab->  Xa,  Xb,  and  Xab  has  been  suppressed  to  simplify  the  notation. 
The  advantage  of  describing  the  voice-call  state  in  terms  of  the  residual  capacity  is  that  the 
number  of  states  is  reduced  to  T  +  1.*'^  Unfortunately,  Cab  is  not  Markovian;  if  it  were,  a 
complete  description  of  the  state  of  the  link  could  then  be  given  by  the  pair  {R,  Cab),  i  e.,  a  pair 
consisting  of  the  number  of  data  packets  in  the  queue  and  the  residual  capacity.  However,  in  our 
approximate  model,  we  make  the  assumption  that  indeed  Cab  is  first-order  Markovian.  Under 
this  assumption,  the  transition  probability  matrix  for  Cab  can  be  obtained  numerically  from  the 
three-dimensional  voice-state  descriptions  =  {xa,  Xb,  Xab),  which  we  approximated  to  be  first- 
order  Markovian  in  our  derivation  of  the  mini-product-form  model  (by  assuming  that  the  three 
flows  are  Poisson  processes).  The  derivation  proceeds  as  follows. 


Each  element  of  the  transition  probability  matrix  for  Cab  can  be  written  as 


Pr(c'  =  C'lc  =  C) 


Pr(c^  =  C^c  =  C) 
Pr(c  =  C) 


where  c  represents  Cab  at  slot  k  and  c'  represents  Cab  at  slot  k+1.  The  denominator,  Pr(c  =  C),  is 
obtained  from  the  equilibrium  distribution  for  s  (i.e.,  Pr(s  =  X))  by  sununing  the  probabilities  of 
each  three-dimensional  voice  state  that  corresponds  to  the  value  c  =  C.  The  equilibrium 
distribution  for  x  is  obtained  by  raising  the  transition  probability  matrix  for  x  to  a  sufficiently 
high  power.  The  numerator  is  obtained  from  the  joint  equilibrium  distribution  of  x  at  two 
consecutive  time  slots,  where 

Pr(x'  =  X',x  =  X)  =  Pr(x  =  X)  Pr(x'  =  X'\x  =  X). 


For  example,  for  T  =  4,  the  number  of  voice  states  is  reduced  from  55  to  5,  and  for  7  =  8  it  is  reduced  from  285  to 
9. 
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Numerical  results,  discussed  in  Section  6,  show  that  this  approximate  model  is  slightly 
optimistic.  At  low  throughput  rates  it  produces  expected  queue  sizes  approximately  2%  lower 
than  that  of  the  mini-product-form  model.  At  high  throughputs  (high  enough  to  produce 
expected  queue  sizes  of  30  packets)  the  discrepancy  is  typically  less  than  10%.  Furthermore,  if 
the  highly  accurate  results  supplied  by  the  three-dimensional  model  are  needed,  the  speed  of  the 
calculation  may  be  increased  by  using  the  one-dimensional  model  to  obtain  an  initial  condition 
for  the  much-slower  three-dimensional  model’s  calculation.  Data-packet  delay  results  based  on 
this  model  are  discussed  in  Section  6.2. 


Fig.  3  —  Example  24-node  network  with  10  source-destination  pairs. 


6.  Performance  results 

We  have  applied  the  mini-product-form  model  to  examples  such  as  the  24-node  network 
shown  in  Fig.  3.  The  shaded  lines  in  the  figure  show  the  7  =  10  voice  circuits  that  we 
considered. ’5  In  this  section  we  discuss  the  accuracy  of  the  mini-product-form  model  in 
characterizing  both  the  voice  process  and  the  data-packet  delay.  In  all  examples,  we  assume  the 
use  of  an  “uncontrolled”  system,  i.e.,  voice  calls  are  admitted  as  long  as  transceivers  are 
available  to  support  them. 


Although  the  mini-product-form  model  has  been  developed  to  improve  the  accuracy  of  data-delay  estimates,  we 
only  need  to  consider  the  voice  portion  of  the  integrated  network  to  evaluate  the  accuracy  of  the  mini-product-form 
model  in  predicting  the  time-varying  behavior  of  the  residual  capacity. 
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6.1  Accuracy  of  the  mini-product-form  model  characterization  of  the  voice  process 

Here  we  evaluate  the  mini-product-form  characterization  of  the  voice  process  on  the 
circled  links  in  Fig.  3.  Clearly,  data-packet  delay  depends  strongly  on  the  residual  capacity  Cab- 
Fox  the  case  of  Ta  =  Tb=  T,  we  consider  the  complement  of  residual  capacity,  namely, 

X  =  T  —  Cab- 

Based  on  Eq.  (2),  we  can  rewrite  this  expression  as 

a:  =  Xab+rxiax{Xa,Xb}, 

and  we  can  interpret  x  as  the  effective  number  of  voice  calls  that  are  currently  being  supported  by 
link  {a,b).  In  the  remainder  of  this  report  we  refer  to  x  as  the  “state”  since  it  is  closely  related  to 
the  voice  state  x  =  [xa,  Xb,  Xab)  defined  earlier.  First,  we  use  the  product-form  solution  to  obtain 
7C(x),  the  equilibrium  distribution  of  x.  We  also  evaluate  the  first  and  second  moments  of  the 
“state- visit  duration,”  E[tix)}  and  E{fi{x)}  respectively,  i.e.,  the  time  spent  in  each  state  x. 
These  quantities  are  evaluated  by  means  of  simulations  of  the  entire  network,  the  mini-product- 
form  network,  and  the  reduced-load  approximation  discussed  in  Section  4.2.  Although  such 
timing  estimates  are  not  used  explicitly  in  our  delay  evaluation,  it  appears  clear  that  an  accurate 
characterization  of  the  time-varying  behavior  of  the  voice  process  is  a  necessary  component  of  a 
good  data-packet  delay  model. 

In  Fig.  4,’^  we  compare  the  exact  state-occupancy  distribution  of  link  (5,13)  of  our 
example  24-node  network  (computed  by  means  of  the  product-form  solution)  to  those  computed 
by  using  the  reduced-load  approximation  and  the  mini-product-form  model  estimate.  Link  (5,13) 
is  centrally  located  and  serves  five  voice  circuits.  We  considered  values  of  p  ranging  from  1  to 
12  (pj  =  p,y  =  1,  10)  for  the  case  of  7  =  8  transceivers  per  node.  In  Fig.  4(a),  the  plot  of 

probability  mass  functions  (pmf  s)  shows  that  the  mini-product-form  model  yields  a  curve  that  is 
very  close  to  the  exact  one;  the  two  curves  can  be  distinguished  only  for  larger  values  of  p  at 
state  X  =  8  (at  which  point  Cs.is  =  0).  The  probability  mass  from  the  reduced-load  approximation 
is  generally  shifted  to  the  right  of  the  exact  result,  thus  predicting  lower  residual  capacity. 


In  our  plots  we  use  continuous  curves  to  represent  discrete  quantities  to  allow  easy  comparison. 
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8  transceivers  per  node 


states  state  jc 


(a)  Distribution  (b)  Weighted  error 

Fig.  4  —  State-occupancy  distribution  on  centrally  located  link  (5,13)  for  various  p  values. 


To  evaluate  the  accuracy  of  the  pmf,  we  consider  the  error  percentage  and  weighted  error 
percentage  as  follows: 


iz(x)  error  percentage  =  100 


7c(j:)-7r(j:)[ 

Kix) 


and 

k(j:)  weighted  error  percentage  =100  %(x)  “  100  |it(x)-7t(A:)|, 

where  n(x)  is  the  estimated  (either  by  the  mini-product-form  model  or  by  the  reduced  load 
approximation)  probability  of  being  in  state  x.  We  limit  our  discussion  of  performance  results  to 
the  case  of  the  weighted  measure  because  relatively  high  percentage  errors  in  infrequently  visited 
states  are  not  expected  to  have  a  significant  impact  on  overall  system  performance.  Figure  4(b) 
shows  that  the  mini-product-form  model  delivers  a  more  accurate  estimate  than  the  reduced-load 
approximation.  With  the  mini-product-form  model,  the  weighted  error  percentage  generally 
increases  with  offered  load  (p).  Apparently  this  is  so  because  contention  for  resources  increases 
with  increasing  load,  thus  causing  blocking  that  cannot  be  anticipated  by  the  mini-product-form 
model.  Although  there  is  not  perfect  agreement,  the  mini-product-form  model  does  provide  a 
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good  estimate  of  the  state  occupancy  distribution  on  link  (5,13);  the  weighted  error  percentage 
shows  that  the  states  where  there  is  a  large  error  are,  in  fact,  infrequently  visited. 

Figure  5  shows  the  voice-call  blocking  probability  on  link  (5,13)  and  the  average  voice- 
call  blocking  probability  for  the  overall  network  (where  the  average  is  taken  over  each  call  type). 
For  the  present  case  of  T  =  8  transceivers  at  each  node,  the  blocking  probability  is  about  0.5  or 
greater  for  p  >  4.  Thus  our  results  for  p  >  4  represent  operation  in  a  regime  of  high  blocking 
probability. 


Fig.  5  —  Blocking  probability  in  the  24-node  network  and  on  link  (5,13); 

8  transceivers  per  node. 

We  continue  our  evaluation  of  the  mini-product-form  model  of  link  (5,13)  in  Fig.  6, 
where  we  plot  the  state-occupancy  distribution  (Fig.  6(a))  and  the  weighted  error  percentage 
(Fig.  6(b))  obtained  from  four  experiments.  The  experiments  differed  in  the  number  of 
transceivers  per  node,  as  shown  in  the  figure.  The  value  of  p  was  fixed  at  one  half  the  number  of 
transceivers  per  node.  For  example,  in  the  experiment  with  four  transceivers  per  node,  we  set 
pj  =  p  =  2,y  =  1,  ...,  10.  In  the  figures,  the  jc  axis  is  the  “fractional  state,”  defined  as  the  real 
state  divided  by  the  number  of  transceivers  per  node.  The  figure  shows  that  the  mini-product- 
form  model  tends  to  yield  more-accurate  results  as  the  number  of  transceivers  per  node  increases, 
although  reasonably  accurate  results  are  achieved  even  for  T  =  2.  This  is  consistent  with 
reduced-load  approximations  that  are  asymptotically  accurate  in  the  limit  of  heavy  traffic  and 
large  link  capacities  (for  a  fixed  ratio  of  traffic/capacity)  [15, 16, 24, 25]. 
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(a)  distribution  (b)  weighted  error 

Fig.  6  —  State-occupancy  distribution  on  centrally  located  link  (5,13)  for  various  numbers  of  transceivers  per  node. 

In  Fig.  7  we  compare  the  first  and  second  moments  (E{t{x)}  and  E{t  {x)}  respectively) 
of  the  state-visit  duration  obtained  from  simulations  of  the  entire  network  (with  eight  transceivers 
per  node,  and  various  p  values),  of  the  mini-product-form  network,  and  of  a  “reduced-load 
network”  (i.e.,  a  two  node,  two  flow  network,  where  the  loads  offered  to  the  two  flows  are  at  the 
rate  given  by  the  reduced-load  approximation).  Here  again  we  see  that  the  mini-product-form 
delivers  more-accurate  estimates  than  the  reduced-load  approximation.*'^ 


To  evaluate  the  accuracy  of  the  mini-product-form  and  reduced-load  models  we  have 
evaluated  mean  error  percentages  and  mean  weighted  error  percentages  for  state-visit  durations 
as  follows: 


and 


mean  £{t(x)}  error  percentage 


1 

T  +  \ 
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£Wx)}-£Wx)}| 

£:{t(x)} 


mean  weighted  error  percentage 


1  |,100  K{x)  E{t(x))-E{t{x)) 


7  +  1 


x=0 


E{t{x)] 


For  values  of  p  of  4  or  greater,  some  of  the  low  occupancy  states  are  infrequently  visited,  thus  resulting  in 
potentially  large  percentage  errors  in  the  simulation  results.  To  reduce  the  possibility  of  artificially  inflating  the 
probability  of  visiting  these  states,  the  initial  state  in  the  simulations  was  two  calls  per  circuit  (rather  than  starting 
from  the  empty  state).  The  simulations  were  run  until  10^  voice-call  arrivals  had  occurred.  In  Fig.  7(a)  the  plot  is 
divided  by  a  stairstep  line  on  7i(a:)  =10"^;  estimates  on  the  left  side  of  this  line  are  based  on  fewer  than  ICP  samples, 
and  are  therefore  subject  to  larger  sampling  error. 
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with  a  similar  definition  for  E{t  (j:)}.  The  quantities  with  the  tildes  are  the  approximate  values 
(based  on  the  mini-product-form  or  reduced  load  approximations),  and  those  without  are  the 
“exact”  values  as  determined  numerically  (7t(j:))  or  by  simulation  (£■{?( jc)}  and  £'{t^(x)}). 


(a)  first  moment  (b)  second  moment 

Fig.  7  —  First  and  second  moments  of  state- visit  duration  on  centrally  located  link  (5,13). 

In  Fig.  8,  we  summarize  the  accuracy  of  the  mini-product-form  model  and  the  reduced- 
load-approximation  estimates  of  the  expected  state-visit  duration  by  plotting  the  mean  error 
percentage  (Fig.  8  (a))  and  the  mean  weighted  error  percentage  (Fig.  8  (b))  as  p  varies  from  1  to 
12.  The  figure  shows  that  the  mini-product-form  model  delivers  significantly  more  accurate 
estimates  than  the  reduced-load  approximation.  Note  that,  in  terms  of  mean  weighted  error 
percentage,  the  mini-product-form  model  is  relatively  insensitive  to  the  network  load,  whereas 
the  accuracy  of  the  reduced-load  approximation  decreases  rapidly  with  increasing  load. 

We  have  examined  in  a  similar  manner  what  we  feel  is  a  representative  set  of  links  in  the 
24-node  network  of  Fig.  3,  and  compiled  the  results  for  easy  comparison  in  Figs.  9  and  10.  The 
links  we  examined  in  addition  to  link  (5,13)  are  link  (15,16),  a  peripheral  link  with  flow/^  =  0, 
three  circuits  deliver  the  traffic  for  flows /ai,  and  fa,  link  (1,2),  a  peripheral  link  that  serves  three 
circuits;  link  (6,8),  a  peripheral  link  that  serves  four  circuits;  link  (1,4),  a  link  that  serves  five 
circuits;  and  link  (5,14),  a  centrally  located  link  that  serves  eight  circuits,  but  has/a*  =  0.  In 
these  figures  we  show  the  mean  weighted  error  percentage  for  the  state-occupancy  distribution 
(Fig.  9),  and  the  first  and  second  moments  of  the  expected  state-visit  duration  (Fig.  10). 
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(a)  mean  error  (b)  mean  weighted  error 

Fig.  8  —  Mean  error  percentage  of  state-visit  duration  on  centrally  located  link  (5,13). 


Fig.  9  —  Mean  weighted  error  percentage  in  state-occupancy  distribution  for  various  links. 

The  results  shown  in  Figs.  9  and  10  suggest  that  the  mini-product-form  model  is  more 
accurate  when  applied  to  links  that  serve  many  voice  circuits.  For  example,  the  mini-product- 
form  characterization  of  the  statistics  for  links  (15,16)  and  (1,2)  is  generally  less  accurate  than 
that  of  the  other  links.  These  support  only  three  circuits  each.  In  addition,  flow/^  on  the  mini- 
product-form  model  of  link  (15,16)  is  zero.  The  results  also  show  that,  as  we  noted  in  our 
detailed  study  of  link  (5,13),  the  accuracy  tends  to  decrease  as  the  offered  load  is  increased. 
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Fig.  10  —  Mean  weighted  error  percentage  in  first  two  moments  of  expected  state- visit  duration 

using  the  mini-product-form  model. 

6.2  Data>packet  delay  evaluation 

Data-packet  delay  has  been  evaluated  using  the  models  discussed  in  this  report.  As 
discussed  in  Section  3,  the  expected  delay  is  obtained  from  the  expected  queue  size  by  means  of 
Little’s  formula.  Figures  1 1  and  12  show  data-packet  delay  as  a  function  of  data-packet  arrival 
rate  for  two  selected  links  of  the  network  of  Fig.  3,  namely  links  (5,13)  and  (1,2),  for  the  case  of 
r  =  4  transceivers  at  each  node  and  for  voice  call  loadings  of  pj  =  =  1  and  2,j  =1,2, ...,  10. 

These  loadings  result  in  overall  blocking  probabilities  (i.e.,  fraction  of  calls  arriving  to  entire 
network  that  are  blocked)  of  0.352  and  0.591,  respectively;  thus  both  cases  represent  a  significant 
amount  of  blocking,  and  hence  a  significant  amount  of  interaction  between  the  links  of  interest 
and  the  rest  of  the  network.  The  blocking  probabilities  for  the  links  of  interest  are 

shown  with  each  figure,  as  are  the  expected  residual  capacities  Cab-  The  expected  voice-call 
duration  (l/p'O  is  100  data  packets  in  all  cases.  The  data-packet  arrival  rate  Xab  is  the 
superposition  of  the  traffic  rate  from  node  a  to  node  b  and  from  node  b  to  node  a.  The 
“utilization”  of  the  residual  capacity  of  link  (a,b)  by  data  traffic  is  then  Xab/Cab^  and  the  delay 
over  link  (a,b)  approaches  infinity  as  Xab  approaches  Cab- 
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a)  p’'  =  2  (C5.,3  =  0.64;  Pb^sm)  =  0.533)  b)  p’'  =  1  (C;.,3  =  1 . 1 7;  PB^5^3)  =  0.3 1 1) 

Fig.  1 1  —  Data-packet  delay  over  link  (5,13)  (7  =  4  transceivers  per  node). 


a)  p''=  2  (C,.2  =  1.67;  Pbo.2)  =  0.138) 

Fig.  12  —  Data-packet  delay  over  link 


data  arrival  rate 


b)  P'^  =  1  (C,.2  =  2.22;  />b(i.2)  =  0.0575) 
,2)  (7  =  4  transceivers  per  node). 


In  Figs.  11  and  12  “exact”  refers  to  results  obtained  from  simulations  of  the  entire 
network,  “mini-product-form”  refers  to  simulations  of  the  three-flow  network,  and  “reduced-load 
approximation”  refers  to  simulations  of  the  two-flow  model  discussed  in  Section  4.2.  All 
simulations  were  run  for  a  duration  of  10’^  data-packet  arrivals.  The  term  “3-dim  Markov” 
(results  currently  available  only  for  link  (5,13)  for  X5  13  <  0.25  and  =  2)  refers  to  the  numerical 
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results  obtained  from  the  Markovian  model  based  on  the  three-flow  mini-product-form  model, 
and  “1-dim  Markov”  refers  to  the  simplified  model  of  Section  5  in  which  the  residual  capacity  is 
assumed  to  be  Markovian.  The  M/D/1  model  is  based  on  a  constant  data-packet  service  rate 
equal  to  the  average  residual  capacity  of  the  link 

Figure  1 1  shows  that  the  mini-product-form  model  accurately  predicts  data-packet  delay 
over  link  (5,13)  at  low-through-moderate  utilization  ranges.  Although  the  accuracy  is  not  as 
good  at  high  utilization,  the  mini-product-form  model  appears  to  provide  a  relatively  close  upper 
bound  on  delay,  which  becomes  looser  as  Xab  approaches  Cab-  Especially  significant  is  the  fact 
that  the  results  for  the  one-dimensional  Markov  model  are  virtually  indistinguishable  from  those 
of  the  simulated  mini-product-form  model,  thus  demonstrating  that  the  reduction  in  complexity 
does  not  come  at  the  expense  of  a  significant  loss  in  accuracy. By  contrast,  the  reduced-load 
approximation  results  in  large  errors,  severely  overestimating  delay  at  moderate  to  high 
utilization  rates. As  expected,  the  M/D/1  model  provides  a  very  poor  (and  extremely 
optimistic)  estimate  of  delay.  The  poor  performance  of  the  Ml/D/1  model  is  easily  explained; 
since  the  expected  voice-call  duration  is  much  greater  than  the  data-packet  length,  the  system  can 
remain  in  a  state  with  low  instantaneous  residual  capacity  for  a  significant  period  of  time,  during 
which  the  data-packet  queue  can  grow  to  large  values;  such  behavior  is  not  predicted  by  the 
M/D/1  model. 

Figure  12  shows  that  the  mini-product-form  model  overestimates  the  delay  for  link  (1,2). 
Although  the  accuracy  of  this  estimate  is  not  as  good  as  that  obtained  for  link  (5,13),  it  is  much 
better  than  that  provided  by  either  the  M/D/1  or  reduced-load  models.  The  fact  that  the  mini¬ 
product-form  model  is  more  accurate  for  link  (5,13)  than  for  link  (1,2)  is  consistent  with  the 
discussion  pertaining  to  Figs.  9  and  10. 

6.2.1  Data-packet  queue-size  distribution 

As  noted  earlier,  our  delay  values  have  been  obtained  by  means  of  Little’s  formula,  which 
relates  the  expected  packet  delay  to  expected  queue  size.  We  have  also  examined  the  distribution 


Although  results  for  the  one-dimensional  Markov  model  are  not  currently  available  for  data-packet  arrival  rates 
greater  than  1.0,  the  simulated  mini-product-form  results  demonstrate  the  accuracy  of  the  mini-product-form 
approach  for  higher  data  rates  as  well. 

These  results  are  representative  in  the  sense  that  the  mini-product-form  model  generally  performs  better  than  the 
reduced-load  model.  Although  the  reduced-load  model  does  provide  smaller  errors  for  some  voice-state  occupancies 
for  some  links  under  particular  loading  considerations  (see  Fig.  4(b)),  in  all  of  the  examples  we  have  studied,  the 
mini-product-form  model  produces  a  more-accurate  estimate  of  data-packet  delay  over  the  entire  range  of  data- 
traffic  arrival  rate. 
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of  queue  size  (rather  than  just  its  first  moment)  in  an  effort  to  understand  better  the  delay 
performance.  Figure  13  shows  the  distribution  of  queue  size^o  for  the  example  corresponding  to 
Fig.  1 1(a),  namely  the  case  of  link  (5,13)  with  =  2  and  an  expected  voice-call  duration  of  100 
times  that  of  data  packets,  for  data-packet  arrival  rates  of  =  0. 1,  0.33,  and  0.5  packets  per  time 
slot.  These  results  have  been  obtained  numerically  using  the  one-dimensional  Markovian  model. 
A  modified  logarithmic  scale  has  been  used  for  the  horizontal  axis  to  illustrate  better  the 
behavior  for  small  values  of  queue  size. 

Figure  13(a)  shows  that  for  =  0.1,  the  queue  is  empty  almost  half  of  the  time  and 
contains  at  most  one  packet  about  60%  of  the  time.  However,  the  probability  mass  at  higher 
values  is  sufficient  to  cause  significant  delay,  as  compared  to  the  results  predicted  by  the  M/D/1 
model.  As  increases,  the  probability  mass  shifts  away  from  the  empty  queue.  However,  it  is 
interesting  to  note  that  the  two  most  likely  queue  sizes  are  still  0  and  1  for  data  rates  at  least  as 
high  as  0.5  (a  relatively  high  data  rate  based  on  the  expected  residual  capacity  of  0.64).  In  all 
three  curves,  there  is  significant  probability  mass  located  at  much  greater  queue  sizes  than  the 
expected  value,  which  results  in  a  large  variance  of  queue  size.  In  fact,  the  standard  deviation  of 
the  queue  length  is  greater  than  its  expected  value,  as  shown  in  Table  1.  The  large  values  of  the 
first  two  moments  of  queue  size  can  be  attributed  principally  to  the  fact  that  the  expected  voice- 
call  duration  is  much  greater  than  the  data  packet  duration,  resulting  in  intervals  during  which  the 
data-packet  queue  size  can  grow  to  large  values,  as  discussed  earlier.  For  this  reason,  the  data- 
packet  array  size  used  in  the  numerical  evaluation  of  the  Markov  chain  must  be  sufficiently  large 
so  that  significant  probability  mass  is  not  lost. 

Table  1  —  First  two  moments  of  queue  size  for  the  example  in  Fig.  1 1(a). 


Xd 

Efqueue  size} 

Std.  Dev. {queue  size} 

0.1 

2.190 

3.386 

0.33 

18.587 

22.970 

0.5 

90.866 

94.171 

20  We  define  the  “queue  size”  at  any  time  to  include  the  packet  being  served  at  that  time. 
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Fig.  13  —  Probability  mass  functions  of  data-packet  queue  size  for  the  example  in  Fig.  1 1(a). 


7.  Conclusions 

A  major  challenge  in  the  development  of  networking  models  is  the  attainment  of  a 
sufficient  degree  of  accuracy  under  the  constraint  of  acceptable  levels  of  complexity.  For 
example,  an  exact  model  for  the  evaluation  of  data-packet  delay  in  integrated  wireless  networks 
(under  the  assumptions  made  in  this  report)  requires  a  7+1  dimensional,  first-order  Markov 
chain,  where  7  is  the  number  of  distinct  voice  circuits.  The  resulting  combinatorial  explosion  in 
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the  size  of  the  state  space  makes  this  exact  approach  infeasible,  except  for  extremely  small 
problems.  On  the  other  hand,  our  simplest  model  for  delay,  the  M/D/1  model,  is  based  on  the 
assumption  that  the  residual  capacity  available  for  voice  is  constant  at  its  expected  value.  Our 
studies  have  demonstrated  that,  not  surprisingly,  such  a  model  predicts  unacceptably  optimistic 
estimates  of  delay  (by  several  orders  of  magnitude  for  typical  system  parameters).  We  have  also 
addressed  the  possible  use  of  a  quasi-static  model,  which  is  based  on  the  fact  that  the  voice  state 
changes  much  more  slowly  than  the  data  state  because  the  expected  voice-call  duration  is 
typically  several  orders  of  magnitude  greater  than  the  data-packet  duration.  However,  we  have 
noted  that  such  models  predict  infinite  delay  whenever  the  offered  data-packet  load  is  greater 
than  the  residual  capacity  in  any  voice  state  with  nonzero  probability.  Thus  such  models  are  of 
no  practical  use  in  estimating  delay  unless  sufficient  capacity  is  reserved  for  the  exclusive  use  of 
data. 


In  this  report  we  have  introduced  a  “mini-product-form”  model  for  the  evaluation  of 
expected  data-packet  delay  in  integrated  multihop  wireless  networks.  In  this  model,  an 
approximate  three-dimensional  Markovian  model  is  used  to  characterize  the  time-varying  voice 
process  at  any  link  of  interest,  resulting  in  significant  reduction  in  complexity  over  the  exact 
system  model.  In  addition,  we  have  developed  a  simplified  one-dimensional  characterization  of 
the  voice  process,  which  very  closely  approximates  the  three-dimensional  model.  Most 
importantly,  we  have  demonstrated,  although  preliminarily,  that  both  of  these  models  can 
provide  accurate  delay  estimates  that  also  happen  to  be  bounds,  at  least  for  certain  networking 
examples.  Further  work  is  needed  to  establish  the  applicability  of  our  models  to  more-general 
scenarios. 

Both  the  mini-product-form  model  and  its  simplified  one-dimensional  version  represent  a 
compromise  between  the  complexity  and  accuracy  of  the  exact  Markovian  model  and  the 
simplicity  but  inaccuracy  of  the  M/D/1  model.  They  provide  better  accuracy  than  the  alternative 
model  that  is  based  on  the  reduced-load  approximation,  which  makes  the  inaccurate  assumption 
that  the  voice  state  at  a  node  is  independent  of  the  voice  state  at  its  neighboring  nodes.  Our 
models  are  certainly  not  the  final  word  on  the  development  of  delay  models.  These  efforts 
represent  a  modest  step  toward  the  development  of  a  universally  applicable  methodology  for 
delay  evaluation.  It  is  hoped  that  they  will  serve  as  the  basis  for  the  development  of  simpler,  but 
acceptably  accurate,  models  that  can  facilitate  the  evaluation  of  system  performance  in  larger 
networks  and  in  more  general  scenarios. 
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